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ABSTRACT 

The gas turbine combustion system design and 
development effort is an engineering exercise to obtain 
an acceptable solution to the conflicting design trade- 
offs between: combustion efficiency, gaseous emis- 

sions, smoke, ignition, restart, lean blowout, burner 
exit temperature quality, structural durability, and 
life cycle cost. For many years, these combustor 
design trade-offs have been carried out with the help 
of fundamental reasoning and extensive component and 
bench testing, backed by empirical and experience 
correlations. 

Recent advances in the capability of computational 
fluid dynamics (CFD) codes have led to their applica- 
tion to complex three-dimensional flows such as those 
in the gas turbine combustor. A number of U.S. Govern- 
ment and industry sponsored programs have made signifi- 
cant contributions to the formulation, development, 
and verification of an analytical combustor design 
methodology which will better define the aerothermal 
loads in a combustor, and be a valuable tool for design 
of future combustion systems. The contributions made 
by NASA Hot Section Technology (HOST) sponsored Aero- 
thermal Modeling and supporting programs are described 
in this paper. 

INTRODUCTION 

The goal of gas turbine combustion system design 
and development is to obtain an acceptable solution to 
the conflicting design trade-offs between combustion 
efficiency, gaseous emissions, smoke, ignition, 
restart, lean blowout, burner exit temperature quality, 
structural durability, and life cycle cost. For many 
years, these combustor design trade-offs have been car- 
ried out with the help of fundamental reasoning and 
extensive component and bench testing, backed by empir- 
ical and experience correlations. The ultimate goal 
has been to develop a reliable combustor design system 
that can provide quantitatively accurate predictions 
of the complex combustion flow field characteristics 
(Fig. 1) so that an optimum combustion system design 
can be achieved within reasonable cost and schedule 
constraints . 


Empirically based procedures have led to success- 
ful evolutionary combustor improvements. However, as 
these methods are experience-based, they are not well 
suited when combustor design requirement are signif- 
icantly different from that of current technology 
engines. The rapidly developing CFD (Computational 
Fluid Dynamics) capability is providing an additional 
tool in the design process which can have a powerful 
positive influence on future design capability. In 
these codes, combustion system subcomponents including 
diffusers, fuel injectors, and combustor liners, in 
addition to the complex internal flow, need to be accu- 
rately modelled. To achieve this, physical sub-models 
and accurate numerical schemes must be developed to 
describe the various aerothermochemi cal processes 
occurring within the combustion chamber. 

A number of U.S. Government and company sponsored 
programs have made significant contributions to the 
formulation, development, and verification of an 
analytical combustor design methodology. These have 
included: U.S. Army Combustor Design Criteria Valida- 
tion (Bruce et al . , 1979; Mongia et al . , 1 979, Mongia 
and Reynolds, 1979), NASA Swirling Recirculating Flow 
(Srinivasan and Mongia, 1980), NASA Soot and NOx Emis- 
sions Prediction (Srivatsa, 1980), NASA Primary Zone 
Study (Sullivan et al . , 1983), NASA Mass and Momentum 
Transfer (Johnson and Bennett, 1981; Roback and 
Johnson, 1983; Johnson et al . , 1984), NASA Lateral Jet 
Injection (Lilley, 1986; Ferrell and Lilley, 1985; 
McMurray and Lilley, 1986; Ong and Lilley, 1986), NASA 
Dilution Jet Mixing (Srinivasan et al . , 1982, 1984, 
1985; Srinivasan and White, 1986; Holdeman et al . , 

1984; Holdeman and Srinivasan, 1986; Holdeman et al . , 
1987a), NASA Transition Mixing Study (Reynolds and 
White, 1986; Holdeman et al . , 1987b), NASA HOST Aero- 
thermal Modeling (Kenworthy et al . , 1983; Sturgess, 
1983; Srinivasan et al., 1983a, 1983b), NASA Error 
Reduction (Syed et al . , 1985), industry IR & D pro- 
grams, and advanced combustor development programs. 

The NASA Hot Section Technology (HOST) Combustion 
Program has supported several of these programs. The 
overall objective of the HOST Combustion Project is to 
develop and verify advanced analytical methods to 
improve the capability to design combustion systems 
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for advanced aircraft gas turbine engines. This objec- 
tive is being approached both computationally and 
experimental 1 y . 

Computationally, HOST first sponsored studies to 
assess and evaluate the capabilities of existing 
aerothermal models (circa 1982). Based on the results 
of these assessments and other studies in the liter- 
ature, HOST supported several studies to develop new 
and improved numerical methods for the analysis of tur- 
bulent viscous recirculating flows, with emphasis on 
accuracy and speed of solution. 

The objectives of HOST sponsored experimental 
studies were to improve understanding of the flow phys- 
ics and chemistry in constituent flows, and to obtain 
ful ly-specif ied, bencbmark-qual i ty experimental data 
suitable for the assessment of the capabilities of 
advanced computational codes. 

This paper reviews the advances in the 
state-of-the-art in combustor aerothermal modeling, 
while highlighting the programs supported by the HOST 
Project (Turbine Engine Hot Section Technology, 1982, 
1983, 1984, 1985, 1986, 1987). Due to length limita- 
tions not all programs that received HOST support are 
included, and, for completeness, some programs that 
made a significant contribution, but which did not 
draw their primary support from HOST are discussed. 

AEROTHERMAL MODELING ASSESSMENT 

Gas turbine combustion models include submodels 
of turbulence, chemical kinetics, turbulence/chemistry 
interaction, spray dynamics, evaporation/combustion, 
radiation, and soot formation and oxidation. A very 
extensive assessment of numerics, physical submodels, 
and the suitability of the available data was made by 
three contractors under Phase 1 of the HOST Aerothermal 
Modeling program (Kenworthy et al . , 1983; Sturgess, 
1983; Srinivasan et al . , 1983a, 1983b). These investi- 
gations surveyed and assessed current models and iden- 
tified model deficiencies through comparison between 
calculated and measured quantities. Results of the 
assessment by Srinivasan et al . , (1983a, 1983b) are 
summarized by Mongia et al . (1986). The constituent 
flows examined included: (1) simple flows with no 

streamline curvature, (2) complex flows without swirl, 
and (3) complex flows with swirl. Geometries for 
several test cases from each of these categories are 
shown in Fig. 2. 

k-e Turbulence Model 

The k-e model is the simplest turbulence model 
that is suitable for recirculating flow calculations. 
This model achieves closure by using a gradient trans- 
port model for Reynolds stress with an isotropic eddy 
viscosity. For flows where the isotropic eddy vis- 
cosity assumption is not valid, the k-e model may be 
either modified (e.g. low Reynolds number correction, 
Richardson number correction) or replaced with an alge- 
braic or differential Reynolds stress model. 

Assessment of the k-e model(s) of turbulence 
showed that these models: 

(1) require low Reynolds Number correction for 
predicting wall shear flows, and streamline 
curvature modifications for accurately 
predicting curved boundary layers 

(2) give quanti tati vely good correlation with 
data for simple flows and non-recirculating 
swi r 1 i ng flows 

(3) give quantitatively reasonable results for 
nonswirling recirculating flows 


(4) give quantitatively unsatisfactory 
correlation with data for complex swirling 
flows with recirculation zones 

(5) give quantitatively unsatisfactory 
correlation, but predict trends correctly, for 
complex three-dimensional flows. 

Algebraic Stress Model and Its Modifications 

Mean flow predictions with this model agreed with 
the data as well as the k-e model results, therefore * 

the conclusions above also apply to this model. In 
addition, the Algebraic Stress Model gives reasonable 
predictions for the Reynolds stress components, con- 
sistent with the strengths and limitations of the k-e * 

models (Mongia et al . , 1986). 

The results of standard k-e and algebraic and dif- 
ferential Reynolds stress turbulence models, have been 
compared in several continuing assessment studies. An 
example comparison (Mongia, 1987) of data and calcula- 
tions using a hybr i d/SIMPLE numerical scheme is shown 
in Fig. 3. This flow is that of co-annular turbulent 
jets flowing into an axisymmetric sudden expansion 
(Roback and Johnson, 1983). In this figure, velocity 
profiles are shown at downstream, distance from 0.11 
to 2.5 pipe diameters from the expansion. 

Scalar Transport Model 

Mongia et al . , (1986) reported that the k-e model 
with specified Prandtl number predicts scalar fluxes 
reasonably well for flow where the gradient diffusion 
approximation is valid. An alternative, the algebraic 
scalar transport model, has the capability to improve 
predictions over the k-e approach, but further work is 
needed to establish its validity for swirling recircu- 
lating flows. 

Turbulence/Chemistry Interaction Models 

It was also concluded by Mongia et al . , (1986) 
that both 2- and 4-step reaction schemes showed prom- 
ise for application in gas turbine combustors, but 
need to be further validated against data from simple 
flames. The modified eddy breakup model predicted 
trends well, and it was recommended that it should be 
pursued because this approach could be easily extended 
to multistep kinetic schemes. 

Numerical Accuracy 

A significant deficiency identified in the assess- 
ments was that for many flows of interest the accuracy 
of the calculation was limited by the numerical approx- 
imations, wherein the false diffusion is of the same 
order of magnitude as the turbulent diffusion. This 
masked the differences between turbulence models such 
that very different models gave essentially the same 
result, and sometimes resulted in undeservedly good 
agreement between data and predictions. 

If false diffusion is present, the numerical solu- 
tion obtained for any given flow depends on the grid 
density and distribution. An example of the compari- 
sons made in the assessment program is given by the « 

comparison in Figs. 4 and 5 between measured and calcu- 
lated temperature distributions downstream from a row 
of jets entering a confined crossflow. This flow is a 
constituent flow in most gas turbine combustors, and i 

has been treated extensively in the literature, includ- 
ing the recently completed NASA Dilution Jet Mixing 
program, from which data were compared with three- 
dimensional calculations in the Phase I assessment 
study by Srinivasan et al . , (1983). 
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The calculated and experimental results shown are 
for a single row of jets with an orifice spacing to 
diameter ratio, S/D, = 2 injected into a ducted main- 
stream with a duct height to orifice diameter ratio 
H/D, = 8. The j et-to-mai nstream momentum flux ratio, 

J, for this test was 25.32. Calculations for this 
case made with 45x26x17 (19890) nodes, are shown in 
Fig. 4. The parameter plotted in these figures is the 
dimensionless mean temperature difference ratio, THETA, 
where THETA = (Tm - T)/(Tm - Tj). The predicted jet 
penetration and mixing are less than that shown by the 
data. 

The calculation shown in Fig. 4 used 49 nodes to 
simulate each jet. It is generally not possible to use 
this many grid points in such a small region; as few 
as four may be used in practice for each jet. To simu- 
late the accuracy of this approximation, calculations 
were performed for the same flow and geometric condi- 
tions, but with a 27x26x8 (5615) grid. These coarse- 
grid calculations (Fig. 5) are in much better agreement 
with the data than the fine-grid calculations. These 
and other calculations in Srinivasan et al . , (1983b) 
clearly demonstrated that the three-dimensional calcu- 
lations were not grid independent. 

Conclusions from the Assessments 

The major conclusion in the HOST Aerothermal 
Modeling Phase I assessment studies by Kenworthy 
et al. (1983), Sturgess (1983), and Srinivasan et al . 
(1983a, 1983b) was that the available computational 
fluid dynamics (CFD) codes provided a useful combustor 
design tool. Although significant advances have been 
made in the development and validation of multidimen- 
sional gas turbine combustion calculation procedures, 
the codes assessed were only qualitatively accurate, 
especially for complex three-dimensional flows, and 
further work was needed. It was concluded that both a 
significantly improved numerical scheme and fully- 
specified experimental data (i.e. both mean and turbu- 
lence flowfield quantities, with measured boundary 
conditions) for complex non-reacting and reacting 
constituent flows were needed before various emerging 
physical sub-models of turbulence, chemistry, sprays, 
turbulence/chemistry interactions, soot formation/ 
oxidation, radiation, and heat transfer could be prop- 
erly assessed. 

A SECOND GENERATION MODEL 

The first generation combustor design procedure 
outlined by Mongia and Smith (1978) has been very use- 
ful for developing several combustors (Mongia et al . , 
1986) that exhibited significant technology advances. 
However, in addition to the model deficiencies identi- 
fied in the assessments, there were several parameters 
of importance in gas turbine combustor design that the 
analytical models could not predict; e.g. gaseous emis- 
sions, soot formation, flame blow-out limits, combus- 
tor pattern factor, and liner heat transfer. These 
parameters were, however, successfully predicted by 
well-established semi -anal yti cal correlations developed 
by Plee and Mellor (1980), Lebfevre (1985), and their 
associates. Therefore, a combustor design procedure 
that could be applied to current and future gas turbine 
engines was implemented that makes use of empirical 
design concepts and employs analytical modeling tools 
to represent various combustion processes (Rizk and 
Mongia, 1986; Mongia, 1987). 

This method makes use of multidimensional models 
to establish liner flowfield features and combustion 
characteristics. The analytical results are then 
integrated with semi -empi ri cal correlations for 
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performance parameters of interest. That is, flow 
field and geometric parameters that are needed in the 
empirical equations, such as combustion volume and the 
fraction of air participating in the primary combustion 
reaction, are provided by the analytical calculations. 

Satisfactory agreement with experimental data has 
been shown (Rizk and Mongia, 1986) for emissions, per- 
formance and heat transfer. The combustor for which 
data were available, and for which calculations were 
performed, is shown schematically in Fig. 6. A typi- 
cal comparison between data and predictions for CO, 
unburned hydrocarbons, NOx , soot emissions, combustion 
efficiency, pattern factor, and lean blowout are shown 
in Figs. 7(a) to (g) respecti vely. The model is in 
good agreement with the data over the entire sea-level 
engine operating range. Calculated liner wall tempera- 
tures for both the inner and outer walls of this com- 
bustor are shown in Fig. 8 for three typical z-planes 
along k = 5, 14, and 23. Here k denotes nodal 
planes along the combustor circumferential direction. 
Although no direct comparison with liner wall tempera- 
ture data was made, the predictions look reasonable. 

AEROTHERMAL MODELING PHASE II 

Based on the recommendations of the Phase I 
assessment studies, activities in Phase II of the HOST 
Aerothermal Modeling program concentrated on developing 
improved numerical schemes, and collecting completely- 
specified data for nonreacting single and two-phase 
swirling and nonswirling flows. The programs initiated 
were: Improved Numerical Methods; Flow Interaction 

Experiment; and Fuel Injector/Air Swirl Characteriza- 
tion. The first of these is a prerequisite to further 
model development, and the data obtained in the latter 
two studies will be used to validate advanced models 
being developed independently. 

Improved Numerical Methods 

The hybrid finite differencing scheme employed in 
generally available combustor codes gives excessive 
numerical diffusion errors which preclude accurate 
quantitative calculations. In response to this defi- 
ciency, HOST supported three programs with the primary 
objective to identify, assess, and implement improved 
solution algorithms applicable to analysis of turbu- 
lent viscous recirculating flows. Both solution accu- 
racy and solution efficiency were addressed (Turbine 
Engine Hot Section Technology, 1985, 1986, 1987; Turan 
and VanDoormal , 1987). 

For most practical problems, a central differ- 
encing scheme would be ideally suited if it were 
unconditionally stable. Central ifferencing is a 
simple second-order scheme which .s easy and straight- 
forward to implement. However, for grid Peclet num- 
bers larger than 2, central differencing can lead to 
over- and under-shoots and is unstable. The hybrid 
(central /upwi nd) scheme is stable for all Peclet num- 
bers, but suffers from excessive false diffusion. An 
alternative scheme, named CONDIF (Controlled Numerical 
Diffusion with Internal Feedback) (Runchal et al . , 

1986) has unconditionally positive coefficients and 
still maintains the essential features of central dif- 
ferencing and its second-order accuracy. 

CONDIF uses central differencing when Pe < 2. 

Where Pe > 2 and the dependent variable varies monoton- 
i cal 1 y , a modified central differencing scheme is 
used, otherwise upwind differencing is used. CONDIF 
employs just enough numerical diffusion to ensure sta- 
bility based internally on the field distribution of 
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the variable, rather than switching to upwind differ- 
encing whenever Pe exceeds 2. Since upwinding is done 
at relatively few grid points, CONDIF essentially main- 
tains the second-order accuracy of central differen- 
cing, and false diffusion is substantially reduced. 

Another advanced numerical scheme, called flux- 
spline (Patankar et al., 1987), is based on a linear 
variation of total flux (convection + diffusion 
between two grid points. This is an improvement over 
the assumption of uniform flux used in hybrid schemes, 
and leads to reduced numerical diffusion. 

Both of these schemes have been used to solve a 
variety of analytical, two-dimensional laminar and tur- 
bulent flows (Runchal et al . , 1987; Patankar et al . , 
1987). As an example, results for a laminar flow 
(Re = 400) in a square driven cavity are shown in 
Fig. 9. This flow, shown schematically in part a), is 
characterized by a strong recirculation zone typical 
of many physical situations. The problem was solved 
with both CONDIF and flux-spine schemes on a uniform 
22x22 grid and compared with the exact analytical solu- 
tion and a hybrid solution on an extremely fine 82x82 
grid. Velocity profiles at the midsection of the cav- 
ity are shown in Fig. 9(b). Both advanced schemes 
show improvement over the hybrid calculation. 

An attractive feature of both CONDIF and flux- 
spline schemes is that their extension to three dimen- 
sions is relatively straight-forward. The resulting 
linear differential equations involve only seven points 
as opposed to 27 points needed in many skewed-upwi nd 
schemes (Syed et al . , 1985). 

In addition to the need for improved numerical 
accuracy, there is a need for improved computational 
efficiency for a given level of accuracy. Typically 
the continuity and momentum equations are solved sepa- 
rately, and then linked through iteration of the 
pressure term; e.g. SIMPLE (Semi -Impl i ci t Method for 
Pressure Linked Equations). Modifications, such as 
SIMPLER and PISO, have been shown to improve computa- 
tional efficiency. Other advanced schemes (Turbine 
Engine Hot Section Technology, 1985, 1986, 1987; Vanka, 
1987), such as block correction techniques and direct 
solution of the coupled equations have been proposed. 
Calculations with the latter coupled with the flux- 
spine technique have shown a speed increase by a fac- 
tor of 15 for a calculation of turbulent flow over a 
backward-facing step (Mongia, 1987). 

Gas Phase Experiments 

An experimental study of the interactions between 
the combustor and diffuser systems (Srinivasan and 
Thorp, 1987) is in progress to: 

(1) Identify the mechanisms and magnitude of 
aerodynamic losses in various sections of an 
annular combustor-diffuser system 

(2) Determine the effects of geometric changes in 
the prediffuser, dome, and shroud on these 
losses 

(3) Obtain a data base to assess current and 
advanced aerodynamic computer models for 
predicting these complex flowfields 

(4) Upgrade the analytical models based on the 
experimental data 

(5) Design and test advanced diffuser systems to 
verify the accuracy of the upgraded analytical 
model 

Another study in progress will obtain comprehen- 
sive mean and turbulence measurements of velocity and 
species concentration in a three-dimensional flow model 
of the primary zone of gas turbine combustion chambers 
(Turbine Engine Hot Section Technology, 1985, 1986, 


1987). The flowfield of interest is the interaction 
between swirling flow and lateral jets in a rectan- 
gular channel (Fig. 10). The mainstreams flow enters 
through 5 swirlers with the transverse jets injected 
from both the top and bottom duct walls with either 2 
or 4 jets per swirler at 1/2 or 1 channel height down- 
stream from the swirler. 

These experiments are being conducted on both air 
and water mul ti pi e-swi r 1 er rigs, as well as single 
swirler and swirling jet rigs. Fifteen cases (combina- 
tions of swirl and jet strength and location) are under 
test using laser sheet light and dye water flow visual- 
ization, and detailed velocity and scalar mean and tur- 
bulence LDV measurements are being made in the air rig. 

A key feature of this program is comparison of 
model calculations against the data obtained to ensure 
that the data are complete and consistent, and satisfy 
the boundary condition input requirements of current 
three-dimensional codes. Calculations were performed 
using a three-dimensional code (Srivasta, 1980) for all 
test cases before the experiments were begun. Data and 
both previous and advanced model calculations are being 
compared as data are obtained. 

Fuel -Injector /Al r-Swi r 1 Characterization 

The objective of this study is to obtain fully- 
specified mean and turbulence measurements of both gas 
and droplet phases downstream of a fuel injector and 
air swirler typical of those used in gas turbine com- 
bustion chambers. 

The flowfield of interest is an axisymmetric 
particle-laden jet flow with and without confinement 
and co-annular swirling air flow. Approximately 30 
cases are under test with both glass-bead particle- 
laden jets and liquid sprays, with various combination 
of swirl strengths and confinement (Turbine Engine Hot 
Section Technology, 1985, 1986, 1987). Measurements 
of mean and turbulence quantities, for both gas and 
solid phases are being made using a 2-component Phase/ 
Doppler LDV particle analyzer (McDonell et al., 1987). 

Calculations were performed for all test cases 
with a two-dimensional TEACH-type nonreacting turbulent 
viscous two-phase flow code before the experiments were 
begun. Data and both previous and advanced model 
calculations are being compared as data are obtained 
(Mostafa et al . , 1987, 1988; Nikjooy et al., 1988). 

In the first series of tests, the developing 
regions of unconfined single and two-phase flows, with 
105 pm glass beads, have been examined experimentally 
and analytically for parti cl e-to-gas mass loadings of 
0.2 and 1.0. Data and calculations for the latter are 
shown in Fig. 11. A two-component Phase/Doppler sys- 
tem was used to map the flowfield, including particle 
number density, and two orthogonal components of veloc- 
i ty for both phases . 

Calculations are shown for both deterministic and 
stochastic treatrents of the particles, using a two- 
phase k-e model. Both treatments of the particles 
give the same gas-phase axial velocity profiles, how- 
ever, the stochastic approach, which attempts to model 
particle/gas phase interactions, gives better agree- 
ment for particle quantities than the deterministic 
approach which ignores turbulence interactions. 

Another experimental program was conducted to 
obtain information on the characteristics of the spray 
produced by a gas turbine fuel injector (McVey et al., 
1988a, 1988b). The objective of this study was to 
obtain spatial ly-resol ved information on both the 
liquid and gaseous phases of the spray flow field under 
conditions of high-flow, high velocity, and high swirl 
that are typical of engine operation. Measurements 
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were made with a hi gh-resol ution spray patternator, a 
two-component laser velocimeter, and a single-component 
Phase/Doppler particle analyzer. 

The comprehensive experimental data generated in 
these programs will be used to validate advanced models 
of turbulence, scalar, and spray transport, including 
two-equation turbulence models, algebraic and differen- 
tial Reynolds stress models, scalar and scalar-velocity 
transport models, and Eulerian and Lagrangian determin- 
istic and stochastic spray models. 

SUMMARY 

Although significant progress has been made in 
the development of three-dimensional analytical CFD 
codes and their application in future gas turbine com- 
bustor design, these codes are neither sufficiently 
comprehensive nor quantitatively accurate enough to 
permit a complete design alone. They are, however, a 
valuable component in an evolving combustor design 
methodology in which their capability is integrated 
with the substantial base of empirical experience and 
one-dimensional flow model ing. 

CONCLUDING REMARKS 

The NASA HOST sponsored Aerothermal Modeling 
Phase II programs will lead to significant improve- 
ments in our technical ability to predict nonreacting 
gas turbine combustor flow fields with and without 
spray injection. Significantly enhanced capabilities 
for accurately predicting combustor aerothermal per- 
formance and wall temperature levels and gradients 
will require further improvements in numerical schemes 
and physical submodels. It is equally important to 
collect fully-specified reacting flow data, similar to 
what is being done for nonreacting flows under HOST 
Phase II, for both complex constituent flows, and 
generic gas turbine combustors. 

In parallel, work should continue in the formula- 
tion and systematic validation of turbulent combustion 
models for reacting sprays and multidimensional heat 
transfer models. These capabilities will provide the 
tools needed to analytically conduct the combustion 
trade-off studies so that optimum future combustion 
systems can be designed, fabricated, and developed 
within acceptable cost and schedule constraints. 
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FIGURE 5. - COMPARISON BETWEEN MEASURED AND CALCULATED DIMENSIONLESS TEMPERATURE DIFFERENCE RATIOS DOWNSTREAM FROM A ROW 
OF COOL JETS INJECTED INTO A CONSTANT-TEMPERATURE CROSS FLOW FROM THE UPPER WALL OF A CONSTANT AREA DUCT (J = 25.32, 
S/H = 0.25, H/D = 8, 5615 NODES). 
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